Transverse Depinning in Strongly Driven Vortex Lattices with Disorder 
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Using numerical simulations we investigate the transverse depinning of moving vortex lattices 
interacting with random disorder. We observe a finite transverse depinning barrier for vortex lattices 
that are driven with high longitudinal drives, when the vortex lattice is defect free and moving in 
correlated ID channels. The transverse barrier is reduced as the longitudinal drive is decreased 
and defects appear in the vortex lattice, and the barrier disappears in the plastic flow regime. At 
the transverse depinning transition, the vortex lattice moves in a staircase pattern with a clear 
transverse narrow-band voltage noise signature. 

PACS numbers: 74.60. Ge 
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The dynamics of driven vortex lattices interacting with 
disorder exhibit a wide variety of interesting nonequi- 
librium behavior and dynamic phase transitions. Ex- 
periments [|l]-||, simulations P^p^, and theory [| 15-1^] 
suggest that at low drives the vortex lattice is disordered 
and exhibits plastic or random flow while at higher drives 
the lattice can undergo a reordering transition and flow 
elastically. In this highly driven state it was suggested 
by Koshelev and Vinokur [p| that the flux lattice forms 
a moving crystal. In subsequent theoretical work, Gi- 
amarchi and Le Doussal lll,|lHl proposed that the re- 
ordered state is actually an ordered moving glass phase 
and that the vortices travel in highly correlated static 
channels. In other work [^jTsI, it has been proposed 
that these channels may be decou pled, p roducing a smec- 
tic structure. Simulations [[lO| , [ll|Jl^Jl4[] and experiments 
|[7| have found evidence for both smectic as well as more 
ordered moving vortex lattice structures. 

A particularly intriguing prediction of the theory of 
Giamarchi and Le Doussal is that, in the highly driven 
phase, the moving lattice has a diverging potential bar- 
rier against a transverse driving force, resulting in the ex- 
istence of a finite transverse critical current. This trans- 
verse critical current has been observed in simulations by 
Moon et al. [Q and Ryu et al. |ll] in the highly driven 
phase. Large transverse barriers have also been seen in 
systems containing periodic pinning |19| . Also, recent ex- 
periments S involving STM images of moving vortices 
reveal that the vortex lattice moves along one of its prin- 
ciple axes rather than in the direction of the drive (as 
predicted by pQl), suggesting that the moving lattice is 
stable to a small transverse force component. 

Although the existence of a transverse critical current 
has been confirmed in simulations, there has been no nu- 
merical study of the properties of the critical current, 
such as the dependence of the barrier size on the strength 
of the longitudinal drive or on the defectiveness of the 
vortex lattice. It would also be very interesting to under- 
stand the dynamics of the vortices at the transverse de- 
pinning transition, and relate this to experimental mea- 



sures such as voltage noise spectra. 

In this work we report a simulation study of the trans- 
verse depinning transition in driven vortex lattices inter- 
acting with random disorder. We find that at high longi- 
tudinal drives, when the vortex lattice is defect free and 
moves in correlated ID channels along one of its princi- 
ple axes, a finite transverse depinning barrier is present. 
For lower drives the transverse barrier is reduced but 
still present, even in the decoupled channel limit when 
adjacent channels slip past each other and some defects 
in the vortex lattice appear. For the lowest drives the 
flux lattice becomes highly defected as it enters the plas- 
tic flow phase, and the transverse barrier is lost. In the 
high driving limit at the transverse depinning transition, 
the vortex lattice spends most of its time moving along 
the longitudinal direction, but periodically jumps in the 
transverse direction by one lattice constant. The vor- 
tex lattice can thus be seen to move in a staircase like 
fashion, keeping its principle axis aligned in the original 
direction of longitudinal driving. As the transverse force 
is increased the frequency of the jumps in the transverse 
direction increases. This motion produces a clear wash- 
board signal in the transverse velocity which can be de- 
tected for transverse drives up to ten times the transverse 
depinning threshold. 

We consider a 2D slice of a system of superconducting 
vortices interacting with a random pinning background. 
The applied magnetic field H = Hi is perpendicular to 
our sample, and we use periodic boundary conditions in 
X and y. The T = overdamped equation of motion for 
a vortex is: 



V^i 



(1) 



where f^ is the total force acting on vortex i, v^ is the 
velocity of vortex z, and rj is the damping coefficient, 
which is set to 1. The repulsive vortex- vortex inter- 



Ar„ 



action is given by f™ = E,=i ^I'/o^idi'i - rjl/A)f 



where r^ is the position of vortex z, A is the penetra- 
tion depth, /o = <I'o/87rA'^, the prefactor Ay is set to 3 
p3|, and Ki{r/\) is a modified Bessel function which 



falls off exponentially for r > A, allowing a cutoff in the 
interactions to be placed at r — 6A for computational 
efficiency. We use a vortex density oi n^ — 0.75/A^ giv- 
ing the number of vortices iV„ = 864 for a sample of 
size 36 A x 36 A. The pinning is modeled as randomly 
placed attractive parabolic traps of radius rp = 0.3A with 

C = (/pAp)(|r. - ri^^ l)e(r, - |r, - 4^) |)f(^\ where ri^) 
is the location of pin k, 8 is the Heaviside step function, 
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The pin density is Up — 1.0/A" and the pinning force 
is fp = l-5/o- The Lorentz force from an applied cur- 
rent J = Jy is modeled as a uniform driving force f^ on 
the vortices in the x-direction. We initialize the vor- 
tex positions by performing simulated annealing with 
fd/fo — 0.0. We then gradually increase fd to its fi- 
nal value by repeatedly increasing fd by 0.004/o and re- 
maining at each drive for 10^ time steps, where dt — 0.02. 
If we increase the drive more rapidly than this, the re- 
ordered vortex lattice that forms at higher drives may 
fail to align its principle axis in the direction of the driv- 
ing. Slow increases in fd always produce an aligned lat- 
tice. Once the final fd value is reached we equilibrate 
the system for an additional 2 x 10^ steps and then be- 
gin applying a force in the transverse direction /^ which 
we increase by O.OOOl/o every lO** time steps. We mon- 
itor the transverse velocities Vy = [l/N^) X)i="i v^ • y to 
identify the transverse critical current. 

In Fig. 1 we show Vy versus the transverse drive /^ 
at longitudinal drives of fd/fo — 1-0 and fd/fo — 3.0 
for a system with a longitudinal depinning threshold of 
Ic / fa ~ 0-5- For fd/ fa — 3.0 the vortex lattice is free 
of defects and the vortices move in well defined ID chan- 
nels as seen from the vortex trajectories in the right in- 
set, in agreement with previous simulations |13|| . In this 
case there is clear evidence for a transverse barrier with 
fc/fc ~ 0.01, approximately 100 times smaller than the 
longitudinal depinning threshold, in agreement with ear- 
lier simulations |10| , [ll| . Thus, the vortex lattice resists 
changing its direction of motion. For fd/fo = 1.0, the 
vortex lattice is highly defected and the ID channel struc- 
ture is lost (as seen in the left inset of Fig. 1). In this 
case the transverse barrier is absent since the lattice has 
no particular alignment and can readily change the di- 
rection of its motion. In the absence of pinning, /^ = 
for all drives fd, as indicated by the top curve in Fig. 1. 

We ran a series of simulations in which the final lon- 
gitudinal drive fd was varied in order to determine the 
dependence of the magnitude of the transverse barrier on 
the magnitude of the longitudinal drive, as well as on the 
density of defects in the vortex lattice. In Fig. 2 we plot 
the resulting transverse depinning thresholds /J' and the 
fraction of six-fold coordinated vortices Pg (calculated 
from the Voronoi or Wigner-Seitz cell construction) as a 
function of fd- For longitudinal drives fd/fo > 1-5, there 
are no defects in the vortex lattice (indicated by the fact 



that Pg « 1.0) and f^ is roughly constant, fy/f^ ^ 0.01. 
Below fd/fo ^ 1-5 defects begin to appear in the vortex 
lattice as adjacent moving channels decouple, and the 
overall vortex lattice develops a moving smectic struc- 
ture. [Of. The transverse critical current /^, which is 
still finite in this phase, becomes progressively reduced 
as more defects are generated. At the lowest drives, 
fd/fo ^ 1-0, the transverse critical force is lost when 
the ID channels are completely destroyed and the vortex 
lattice enters the amorphous plastic fiow phase shown in 
the left inset of Fig. 1. The dislocations in the lattice, 
which were aligned perpendicular to the vortex motion 
at fd/fo > 1-0, become randomly aligned at the transi- 
tion to plastic flow. The loss of f^ thus coincides with 
the loss of alignment of the defects and the destruction 
of the ID channels. 

We have also checked the effect of finite system size on 
the magnitude of the transverse critical force for fd/fo = 
1.5,2.0 and 3.0 for different system sizes {L — 24A, 36A, 
48A and 60A). In the inset of Fig. 2, we show that fy / f^ 
is not affected by the system size. These results support 
the idea that it is the presence of defects in the lattice 
that reduce or destroy the transverse barrier, rather than 
any type of matching effect with the system size, and that 
as long as some form of channeling occurs the barrier will 
still be present. 

In order to view the dynamics of the vortex lattice 
at the transverse depinning threshold we plot in Fig. 3 
the vortex positions and trajectories for the same system 
shown in Fig. 1 with fd/fo = 3.0 and f^/f^ = 0.011. 
For clarity we have highlighted a particular row of vor- 
tices. In Fig. 3(a) the principle axis of the vortex lattice 
is aligned with the direction of the drive and the ordered 
lattice is moving along this axis. In Fig. 3(b) the en- 
tire vortex lattice has translated by one lattice constant 
in the transverse direction. During the transition the 
lattice moves at an angle to the longitudinal drive, fol- 
lowing a different axis of the lattice. Once the vortices 
have moved one lattice constant transverse to the drive, 
they begin moving along the same channels that formed 
before the transverse translation. After this, the vor- 
tices move along the longitudinal direction once again, as 
in Fig. 3(c), before jumping by another lattice constant 
in the transverse direction. At the transverse depinning 
transition, the vortex lattice thus moves in a staircase 
like manner, always keeping its principle axis aligned in 
the direction of the original longitudinal drive and always 
translating along one of the axes of the lattice. As /^ is 
increased the frequency of jumps in the transverse direc- 
tion increases. If /J' is increased to a high enough value, 
we have found evidence that the vortex lattice will reori- 
ent itself with the net driving force via the creation of 
a grain boundary. This will be discussed in more detail 
elsewhere. The transverse depinning transition is unlike 
the longitudinal depinning transition in that the latter 
occurs through plastic deformations of the lattice and 



the generation of a large number of defects. In contrast, 
the transverse depinning transition is elastic. 

The fact that the vortices move periodically by a lattice 
constant in the transverse direction is a result of the fact 
that the longitudinal channels followed by the vortices 
are uniquely determined by the underlying disorder [^ . 
The vortices jump from one of these stable channels to 
another, giving the same effect as a washboard potential. 
This periodic effect occurs only for a moving lattice in 
which ID channels have formed; to a stationary lattice, 
the disorder would appear random. 

A consequence of the staircase-like vortex motion just 
above the transverse depinning threshold is that the net 
transverse vortex velocity at a fixed /^ should show a 
clear washboard frequency which should increase for in- 
creasing /|. In Fig. 4(a) we plot Vy for samples with /| 
held fixed at several different values just above the trans- 
verse depinning threshold for a system with fd/fo = 3.0. 
For f^/f^ = 0.011 the Vy shows periodic pulses which 
correspond to the correlated transverse jumps of the vor- 
tex lattice seen in Fig. 3. The fiat portions of the voltage 
signal correspond to time periods when the lattice is mov- 
ing only in the longitudinal direction, between hops. For 
increasing transverse drive the frequency of these pulses 
also increases. The additional structure in the Vy volt- 
age pulses at lower values of /^ is characteristic of the 
underlying pinning, and varies for different disorder real- 
izations. It occurs when the vortex lattice moves slightly 
unevenly, with a small wobble, but no defects or tearing 
occur in the lattice. The main feature of large periodic 
pulses is always observed, and the wobble dies away at 
larger transverse drives. In Fig. 4(b) we show that the 
Fourier transform of the velocity signal Vy for a driving 
force oi fj/fc — 0.016 exhibits a resonance frequency at 
ly = 6.0 X 10^^ inverse MD steps. In Fig. 4(c), the reso- 
nant frequency increases linearly with f^. We find that 
this resonance persists for /^ up to ten times larger then 
the transverse depinning threshold. It should be possi- 
ble to detect this washboard frequency with Hall-noise 
measurements. 

In recent experiments employing an STM to directly 
image a slowly moving vortex lattice ||21| , evidence for 
staircase-like motion of the flux lattice has been observed. 
In these experiments the direction of the driving force 
could not be directly controlled, but was assumed to be 
at a slight angle with respect to the principle vortex lat- 
tice vector, so that a transverse component of the driv- 
ing force was present. Further experiments in which the 
magnitude and direction of the drive can be directly con- 
trolled are needed; however, experimental imaging tech- 
niques such as STM or Lorentz microscopy 123] seem 
highly promising. 

In summary we have investigated the transverse depin- 
ning of moving vortex lattices interacting with random 
disorder. We find that for high longitudinal drives where 
the vortex lattice is defect free a finite transverse bar- 



rier forms. For lower drives where defects in the vortex 
lattice form and the vortex lattice has a smectic struc- 
ture the transverse barrier is reduced but still finite. In 
the highly disordered plastic fiow phase the transverse 
barrier is absent. The transverse depinning transition is 
elastic, unlike the plastic longitudinal depinning transi- 
tion, and near this transition the vortex lattice moves in a 
staircase-like fashion. We observe a washboard frequency 
in the transverse voltage signal which can be detected for 
transverse drives up to ten times the depinning drive. 
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FIG. 1. Transverse voltage Vy versus transverse driv- 
ing force /^. Diamonds: A sample with no pinning, 
fp/ fo ~ 0, shows no transverse barrier at any drive (here, 
fd/fo ~ 3.0). Squares and circles: a sample with pin- 
ning strength fp/fo = 1.5 and longitudinal critical current 
fc/fo = 0.5. Squares: For a longitudinal driving force 
fd/fo = 1.0 the vortices flow plastically (left inset) and there 
is no transverse barrier. Circles: For a longitudinal driving 
force fd/fo = 3.0 the vortices flow in well-deflned channels 
(right inset) and a transverse barrier of fc/fc ~ 0.01 ap- 
pears. Insets: Filled circles represent vortices and lines show 
the paths traveled by the vortices while moving over randomly 
spaced pinning sites (not shown). Left inset: plastic flow; 
right inset: channel flow. 
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FIG. 2. Circles: Transverse critical force fc/fc versus lon- 
gitudinal driving force /c;. Diamonds: Corresponding fraction 
of six-fold coordinated vortices, Pq. The transverse critical 
force saturates for a reordered vortex lattice, fd/fo > 1-8, 
but drops as the lattice becomes defected, reaching zero at 
fd/fo ~ 1.0. Inset: Transverse critical force for different sys- 
tem sizes L = 24A, 36A, 48A, and 60A. Circles: fd/fo = 1-5; 
Squares: fd/fo = 2.0; Diamonds: fd/fo = 3.0. 




(a) (b) (c) 

FIG. 3. Consecutive simulation images of vortex motion 
just above the transverse critical force, f^/fc = 0.011. Cir- 
cles represent vortices and lines indicate paths followed by the 
vortices. A particular row of vortices has been highlighted. 
In (a) the vortices are moving in the direction of the applied 
longitudinal drive. In (b) the vortex lattice changes its di- 
rection of motion and follows a different lattice vector until 
it has translated by one lattice constant. In (c) the lattice 
again switches direction and continues to flow in the longi- 
tudinal channels. The frequency at which the vortex lattice 
hops from channel to channel increases with increasing trans- 
verse drive. 
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FIG. 4. (a) Transverse voltage signal Vy as a function 
of time for a sample with longitudinal drive fd/fo = 3.0 
at four different transverse driving currents /^. Bottom to 
top: f^Jf^ = 0.011, f^Jf^ = 0.014, f^Jf^ = 0.016, and 
fd/fc ~ 0.20. Each pulse corresponds to the vortex lattice 
moving one lattice constant in the transverse direction, (b) 
Power spectrum S{i') of the transverse voltage noise signal 
Vy for a sample with fa/fo = 3.0 and f^/f^ = 0.016. A 
clear narrow band signature appears ai v — 6.0 x 10^^. (c) 
Location of the narrow band peak f for different transverse 
driving forces /^. 



